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1.0  INTRODUCTION 


In  the  contract  proposal  the  Principal  Investiga- 
tors set  out  objectives  for  their  research  which  involved 
finding  new  techniques  for  solving  problems  in  radiative 
transfer  and  neutron  transport.  We  also  proposed  to  start  an 
investigation  of  the  inverse  problem  which  has  only  been 
taken  up  in  recent  years.  The  direct  problem,  which  is  the 
usual  problem  treated,  attempts  to  solve  the  transport  equa- 
tions when  the  phase  function  (which  characterizes  the  emis- 
sion and  absorption  properties  of  the  atoms  or  molecules  in 
the  medium)  is  known.  By  contrast,  the  inverse  problem  is 
concerned  with  finding  the  phase  function  when  sufficient  in- 
formation is  given  about  the  gross  transmission  and  absorp- 
tion properties  of  the  medium  as  a whole.  As  indicated  in  the 
proposal,  our  initial  approach  has  been  directed  to  one-dimen- 
sional problems,  i.e.  problems  in  which  stratification  occurs 
along  one  of  the  axes,  say  the  z-axis.  We  are  happy  to  re- 
port that  we  have  made  considerable  progress  in  our 
objectives. 


2.0  RESULTS  OF  RESEARCH 
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The  results  of  our  research  is  rather  fully  docu- 
mented in  three  papers,  two  of  which  have  been  published  and 
the  third  of  which  is  in  the  process  of  publication.  We  be- 
lieve that  our  achievements  are  considerable,  considering 
that  they  were  completed  within  one  year.  The  inverse  prob- 
lem, particularly,  has  broken  new  ground  and  offers  the  possi- 
bility of  characterizing  atoms  or  molecules  of  a medium  when 
the  transmittance  or  albedo  are  known.  Moreover,  it  also 
suggests  the  possibility  of  a rational  method  of  designing 
nuclear  piles  corresponding  to  neutron  flux  requirements. 


^ . 


D 

0 

0 

0 


u 

LI 


3.0  PAPERS  FOR  PUBLICATION 


As  mentioned  above,  the  accomplishments  are  given 
in  three  papers,  two  of  which  are  published  and  the  third  of 
which  has  been  accepted  for  publication.  Two  reprints  cor- 
responding to  the  published  papers  and  a preprint  of  the 
paper  to  be  published  are  given  as  part  of  this  final  report. 
In  the  present  section  we  give  only  the  abstracts. 


(a)  "A  Variational  Principle  for  Transport  Theory,”  Journal 
of  Mathematical  Physics,  1£,  1258  (1978). 

ABSTRACT 

A maximal  variational  principle  is  used  to  construct  an 
infinite  medium  Green's  function  for  treating  the  bound- 
ary value  problems  of  the  linear  transport  theory  (neu- 
tron and  radiative).  For  the  neutrons  we  consider  the 
one-speed  case  and  correspondingly  for  the  radiative 
transfer  the  monochromatic  case.  The  scattering  proper- 
ties of  the  medum  are  presumed  to  be  dependent  on  the 
relaxation  length.  Thus,  for  the  neutrons  the  secondary 
production  function  depends  on  the  neutron's  relaxation 
length  and  for  the  radiative  transfer  the  albedo  for 
single  scattering  is  dependent  on  the  optical  depth. 


These  two  functions  are  kept  arbitrary  so  that  a large 
class  of  problems  can  be  covered.  The  basic  principle 
involves  a functional  which  is  an  absolute  maximum  when 
the  trial  function  is  an  exact  solution  of  an  integral 
equation  of  the  Fredholm  type.  The  kernel  of  the  inte- 
gral equation  is  required  to  satisfy  certain  symmetry 
and  boundlessness  properties.  We  also  exhibit  an  inter' 
esting  relation  between  the  absolute  maximal  and 
Schwinger's  stationary  variational  principles,  which  in 
general  is  neither  a maximum  nor  a minimum. 


(b)  "Direct-Inverse  Problems  in  Transport  Theory.  1.  The 
Inverse  Problem,"  Journal  of  Mathematical  Physics,  19, 
1793  (1978). 

ABSTRACT 

For  the  inverse  problem  treated  here  we  use  the  results 
of  an  experiment  which  measures  the  total  angular-de- 
pendent column  density  as  compared  to  the  measurement 
which  provides  information  on  the  angle  integrated 
spatial-dependent  angular  density  (or  the  specific  in- 
tensity). We  use  two  methods  of  approach.  One,  the 
Legendre  expansion  method  and  two  a maximal  variational 
principle.  In  particular  we  demonstrate  how  the  varia- 
tional principle  yields  a very  convenient  representation 
of  the  scatterirg  kernel  (or  the  phase  function)  in 
terms  of  a basis  consisting  of  Case  eigenfunctions  for 
the  isotropically  scattering  medium. 


(c)  "Direct-Inverse  Problems  in  Transport  Theory.  The  In- 
verse Albedo  Problem  for  a Finite  Medium,"  Journal  of 
Mathematical  Physics  (to  be  published). 

ABSTRACT 

In  this  paper  we  deal  with  the  inverse  problem  in  radia- 
tive transfer,  which  is  equally  applicable  to  the  neutron 
transport,  for  a finite  homogeneous  medium.  We  give  a 
method  for  computing  the  scattering  function  and  the 
albedo  for  single  scattering.  The  solution  is  given  in 
terms  of  Legendre  expansion  of  the  scattering  function. 

A decomposition  of  the  equation  of  transfer  is  also  given 
in  which  the  relation  between  the  direct  and  the  inverse 
problem  is  exhibited  via  the  principles  of  invariance. 

A relation  of  this  work  with  Case's  method  of  approach  is 
outlined.  This  work  should  be  of  practical  value  to  the 
problems  associated  with  remote  sensing  of  the  terres- 
trial atmosphere  and  the  neutron  piles. 
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For  the  invene  problem  treated  here  we  use  the  results  of  an  experiment  which  measures  the  total  angular- 
dependent  column  density  as  compared  to  the  measurement  which  provides  information  on  the  angle 
integrated  spatial-dependent  angular  density  (or  the  specific  intensity).  We  use  two  methods  of  approach. 

One,  the  Legendre  expansion  method  and  two  a maximal  variational  principle.  In  particular  we 
demonstrate  how  the  variational  principle  yields  a very  convenient  representation  of  the  scattering  kernel 
(or  the  phase  function)  in  terms  of  a basis  consisting  of  Case  eigenfunctions  for  the  isotropically  scattering 
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where,  c < 1 ( nonmultiplying  medium ),  ♦ is  the  distribution 
function  and  Alt' -*p)  is  the  normalized  symmetric  scatter- 
ing function,  i-e.. 


1.  INTRODUCTION 

The  object  of  the  direct  problem  of  transport  theory, 
(eg.,  neutron  or  radiative)  the  usual  problem  treated,  is  to 
find  the  distribution  function  for  a given  scattering  function. 
By  contrast,  the  inverse  problem  involves  construction  of  the 
scattering  function  from  the  results  of  some  simple  experi- 
ment which  gives  some  knowledge  of  the  distribution  func- 
tion. For  the  inverse  problem  in  general,  it  is  usually  not 
clear  with  regard  to  the  nature  of  the  minimal  set  of  measure- 
ments one  needs  in  order  to  construct  the  scattering  function 
uniquely.  For  instance,  consider  a relatively  simple  case  of 
transport  of  monoenergetic  neutrons  in  a medium  with  a 
homogeneous  mixture  of  different  nuclear  species.  If  the  an- 
gular scattering  differential  cross  section  associated  with 
each  species  is  different,  then  it  would  appear  to  be  difficult 
to  devise  an  experiment  from  which  any  set  of  measurements 
can  lead  to  a unique  conclusion  for  the  shape  of  the  angular 
differential  cross  section  for  each  nuclear  species.  If  one  had 
a medium  consisting  of  only  one  nuclear  species,  which  is  not 
necessarily  homogeneously  distributed,  then  certain  rela- 
tively simple  measurements  will  lead  to  a unique  determina- 
tion of  the  differential  cross  section.  The  reason  is  simply 
that  in  that  case  the  secondary  production  function,  which  is 
the  ratio  of  the  rate  at  which  the  neutrons  are  elastically 
scattered  to  the  rate  at  which  the  scattering  involves  all  nu- 
clear processes,  become  independent  of  the  density  of  that 
nuclear  species  and  depends  only  on  its  scattering  properties. 
Similar  considerations  apply  to  problems  of  radiative  trans- 
fer and,  inter  alia,  the  problem  of  transport  of  hyperthermal 
electrons  in  a plasma  containing  a substantial  amount  of 
neutral  species  of  various  kinds. 

2.  DIRECT  INVERSE  PROBLEMS 

For  our  presentation  of  the  direct  inverse  problems  in 
transport  theory  we  consider  an  infinite  medium  with  a 
plane  source  at  x = x.  emitting  ooe-’peed  neutrons  in  a direc- 
tion whose  cosine  is  p*.  For  the  onenlimensional  time-inde- 
pendent, and  azimuthally  symmetric  case,  the  standard  neu- 
tron transport  equation  is1 


ijV/V-lO-l. 


For  the  direct  problem,/!)!’— >|i)  is  given  and  ¥ is  to  be  deter- 
mined everywhere.  For  the  inverse  problem  certain  results  of 
an  experiment  are  given  and  then is  to  be  construct- 
ed from  those  measurements. 

In  an  earlier  paper,*  Case  presents  an  elegant  method  of 
solving  the  inverse  problem.  The  essence  of  his  procedure  is 
that,  if  one  expands  the  scattering  function  in  terms  of  Le- 
gendre polynomials  so  that 


Am'-m)«§(2/+ 1 yrwpfyL). 


then  for  the  unit  isotropic  plane  source  the  expansion 
coefficients 

*<0=1 -c/,  (4) 

are  determined  from  the  measurement  of  the  density  of  neu- 
trons as  a function  of  x.  Therefore,  as  Case  has  shown  from 
the  direct  problem,  one  has  the  spectral  representation  of  the 
total  density  [$(x)J 


*{x)=fc^re ~Wv’ 
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where  pi  v)  is  the  spectral  function  given  by 

= -1<V<1 

V AT(v) 

■C,  8(v-v()</v  , , . 

= 1 jy- . M>1.  (6) 

where  N(v)  and  N,  are  the  normalizations  of  the  eigenfunc- 
tions of  the  continuous  and  discrete  spectra,  respectively. 
For  the  inverse  problem  he  states  that,  given  the  measure- 
ment of  the  total  density  4>(x),  one  knows  quite  a bit  about 
the  spectral  density  from  Eq.  ( 5 ).  We  will  not  give  further 
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details  here  as  to  how  this  argument  actually  leads  to  the 
determination  of  the  expansion  coefficients  g</)  and  hence 
(see  Refs.  2 and  3 for  a complete  discussion). 

We  wish  to  present  a different  method  in  the  present 
paper  for  dealing  with  the  inverse  problem.  In  contrast  to  the 
data  used  in  Case’s  approach  we  require  the  measurement  of 
the  cngular  column  density  M0(p,pe)  of  neutrons  for  a unit 
plane  source  emitting  neutrons  in  one  direction  p,  (mono- 
directional  ).  In  other  words,  the  experiment  should  provide 
the  information  on  the  quantity 

AS«0i,p,)=*J  dx  ¥(x,p-.p,)  (7) 

where  ¥(x,p—.p«)  is  the  angular  density  of  neutrons  satisfy- 
ing Eq.  ( 1 ). 

Clearly  MJj. t,p,)  is  the  angle-dependent  zeroth  mo- 
ment of  the  distribution  function  (the  angular  density) 
which  is  related  to  the  scattering  function.  That  relation  is 
readily  obtained  from  Eq.  ( 1 ) by  integration  with  respect  to 
x.  With  the  appropriate  boundary  conditions  that  ¥ vanish- 
es at x**±ao,  we  have  from  Eq.  (1) 

MfaH,) = 6(p-p«)  + yj 

(8) 

This  is  a Fredholm  type  of  an  integral  equation  for  3/o(p,p,) 
which  involves  the  scattering  function  as  the  kernel.  If  we  let 

5«0i,p,)=J  dp'./(p'-*p)A/,(p',p,)  (9) 

then  S» (fi>M«)  ( the  emission  term ) also  satisfies  a similar  type 
of  an  integral  equation.  That  equation  is  readily  obtained  by 
multiplying  Eq.  ( 8 ) by/0i-»p,)  and  integrating  with  respect 
to  p.  With  change  of  names  of  variables  we  have 

W-"P)+  -jf  dp'Ap'—WW.Pa). 

(10) 


For  convenience  we  write  Eq.  (8)  as 


(11) 

Before  addressing  the  Direct  Inverse  problems,  we  wish  to 
point  out  that  if  in  the  integral  equations  (8)  and  ( 10)  the 
kernel /T p-*p»)  satisfies  certain  properties  then  there  is  a 
maximal  variational  principle*  which  can  be  used  to  solve 
such  equations.  Explicitly,  if/Tp-»p«)  satisfies  the  following 
properties: 

(a) /Tp-.po)  is  symmetric  in  p,p* 

(b) /lp-»p«)>0,  for  — l<(p,p«)<l, 

(c)  -—-  f ji.)<l,  -l<p,<l,c<l, 

2 J- 1 

then  the  functional 


J *.»*•) 

+,y 


(12) 


is  an  absolute  maximum  if  and  only  if  />(p,p«)  is  an  exact 
solution  of  the  integral  equation  ( 10).  Or  conversely,  as  in 
the  inverse  problem,  if  we  consider  Eq.  ( 10)  as  an  integral 
equation  for /X  p«-»p)  so  that 

Ap*~*p)  =S,(p,p,)-  yj  fy'SaHn'.Ha)/ (p'-.p), 

(13a) 

then  the  functional 

FJin)-  J dp  «(p»,p)|  2Sdji,Ha)-n(jia,ft) 

- y £ dp'  So(p',p,)n(p'.p)]  ( 13b) 

is  an  absolute  maximum  if  and  only  if  n ( p*p ) =fl  p«— »p ) is 
the  exact  solution  of  the  integral  equation  (13a).  This  result 
follows  from  the  fact  that  5,(p,pe)  is  symmetric  and  the 
operator  A corresponding  to  the  integral  equation  ( 1 3a ) , i.e. 

\f=f  (13c) 

is  positive  definite.  In  two  previous  papers  by  Kanal  and 
Moses’-*  such  a variational  principle  was  used  to  solve  the 
problems  of  inverse  scattering  and  also  a demonstration  of 
its  application  to  the  linear  transport  theory  for  the  isotropic 
scattering  inhomogeneous  media  was  given.  Here  for  the  di- 
rect inverse  problems  we  can  proceed  in  two  ways.  One  is  to 
expand/Tp— »p«)  in  terms  of  Legendre  polynomials  and  use 
their  orthogonality  property  to  relate  the  expansion  coeffi- 
cients to  the  Legendre  moments  of  So(p>p«)-  This  method 
produces  exact  solutions.  The  other  way  is  to  use  the  vari- 
ational principle.  Use  of  either  method  will  be  dictated  by  the 
kind  of  an  experiment  that  is  required  to  relate  the  scattering 
function /(p-»po)  to  the  emission  term  .SUp.Po)  or  vice 
versa.  We  shall  illustrate  both  methods. 


A.  Inverse  Problem  by  Legendre  Expansion 

First  we  note  that  if/(p— .pj  is  symmetric  then  so 
is  5e(p,p«).  Let  us  now  expand/Tpr-»p«)  in  terms  of  Le- 
gendre polynomials  so  that 


Ap-+p,)  = ^(2,+iitfM.W. 

Insert  this  expansion  in  Eq.  ( 10)  to  obtain 
S,(vi,na)=£  (2n+l)/^(p)F„(M,). 

M-0 

+ 7 £ (2n  + l)/SM)lM 

+ If— 0 

where 

9.0*«)=J  dp  P^Wp,p*)- 


(14) 


(15) 

(16) 


Using  the  orthogonality  property  of  Legendre  polynomials 


we  conclude  from  Eq.  ( 1 5 ) that 
, 1<P q) 

' IPfoiaj+CqtHa) 


(17) 


(18) 
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Note  that  /are  independent  of  po.  In  consequence  one  needs 
to  make  the  measurement  for  any  one  direction  p„  to  deter- 
mine all  the  expansion  coefficients  In  particular,  for 
Po=  1,  i.e.,  for  the  source  emitting  neutrons  parallel  to  the  x 
axis,  we  have 

, 9/0) 

Jl- 


2+cq{l) 


(19) 


This  result  is  really  not  surprising  for  a homo^e-  *ous  medi- 
um due  to  the  fact  that  the  medium  is  rotationalr,  invariant. 
However,  it  is  of  prime  importance  for  the  construction  of 
/[p-*pc).  In  particular,  if  the  medium  is  not  drastically 
anisotropically  scattering,  then  one  needs  only  a few  expan- 
sion coefficients  and  this  method  becomes  very  useful. 

Now  the  Legendre  moments  q, (p,)  ofS0(p,p<,)  can  be 
determined  from  the  measurement  of  the  angular  column 
density  Mj(p,pJ.  In  other  words,  given  MJp.p,),  we  con- 
clude from  Eq.  ( 1 1 ),  that 


9/0*«)=  y [ J , dp  P{p)M°(p,P<>) - ^/Wj • ( 20) 


The  quantity  c ( the  secondary  production  term ) is  obtained 
from  Eq.  (8)  [or  from  Eq.  (11)]  by  integration  with  respect 
to  p so  that 


:=1  — [J  dpM^p.p^) 


(21) 


This  is  a well-known  result  (cf  Ref.  1 ).  However,  from  the 
inverse  problem  point  of  view  we  again  note  that  c,  being 
independent  of  p*,  is  determined  from  the  experiment  with 
the  source  emitting  neutrons  in  any  arbitrary  direction. 

It  would  be  interesting  to  know  the  degree  of  anisotropy 
of  the  scattering  medium.  We  can  estimate  that  from  the 
experiment  if  q,(p<>)  were  calculated  for  all  fi*.  i.e.,  if  we 
calculate  all  the  Legendre  moments  of  the  column  density 
[see  for  example  Eq.  ( 19)].  For  in  that  case  we  obtain  from 
Eq.  ( 18)  the  relation 

dpPjtp)=* 


I 


SnJ. 


(22) 


-i 


(2n+lKl-c/n) 

In  other  words,  (P„(p ),  q/p))  form  a biorthogonal  set,  but 
more  importantly  for  l=n  we  get 

(2* +D7; 


fn 


4+c(2n+ 1)7; 


(23) 


where 


7>J'  d^Oifc.O*). 


(24) 


The  f„'i  will  give  a measure  of  the  degree  of  anisotropy  so 
that  /„= 0,  for  n>N. 


B.  Inverse  Problem  by  the  Variational  Principle 


Expansion  in  terms  of  Legendre  polynomials  is  useful 
when  the  situation  is  not  too  anisotropic.  When  one  wishes 
to  consider  very  anisotropic  cases,  the  variational  principle 
may  be  more  useful.  Thus,  consider  Eq.  ( 13a)  so  that 
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/lp-*Po)=S,(p.p,)-  f dp' S»(p’,p,)/(p'-*p). 
2 J- 1 


(25) 


l{n(p,pt)=f(p-*pi>)  is  the  exact  solution  of  Eq.  (25),  then 
the  functional  defined  by  ( 13b)  is  an  absolute  maximum. 
The  value  of  the  functional  is 


Fj\f]  = — [S0(^lo,/io)-/(jl<r->Mo)] 
c 


(26a) 


or 


/N(pc-<-pt)^So(po,Po)-  y />„[/■*]. 


(26b) 


where  fN  is  a class  of  trial  functions.  From  ( 26a ) we  see  that 
from  the  class  of  functions  {/*},  the  function  which  maxi- 
mizes the  functional  ( 13b)  is  the  desired  one.  From  (26b) 
one  obtains  an  upper  bound  for/t ft*— »ji*),  i.e.,  forward  scat- 
tering for  any  fN  used.’ 

To  illustrate  the  application  of  the  variational  principle 
we  now  give  an  example.  In  analogy  with  the  comparison 
potential  technique  in  the  inverse  scattering  theory,'  we  shall 
choose  a sequence  of  trial  functions  forf(p-+p*)  which  is 
constructed  from  the  complete  set  of  eigenfunctions  of  the 
transport  equation  with  isotropic  scattering.  In  this  way  some 
very  elegant  results  are  obtained.  Thus,  we  .nail  expand 
/lp-*Po)  in  terms  of  those  eigenfunctions  and  find  that  the 
expansion  coefficients  satisfy  decoupled  integral  equations 
analogous  to  Eq.  ( 25 ) involving  the  same  kernel  S0( p.p<>). 

For/Tp.-*/i*)  = 1,  we  have  the  completeness  relation, 


1 


5(ji-ji*)=  — [p.4»*.(f i)<Mfi*)-^<Mjt)<MM«)] 
A'©* 


dv 


N(V) 


(27) 


where  <t>0± (p),  are  Case’s1  discrete  and  continuum 
eigenfunctions,  respectively,  defined  by 


<^o±(M')==± 


CVo 


I 


4>Jn)=jP 


2 ±v,-ft 
1 


(28a) 


v-p 

\(v)  = i[A’(v)  + A(v)], 


+ X(v)8(v— p), 


=1- — r* 

2 z-p 

>)  =0, 


A(Z) 
A(i/„)=0, 


(28b) 

(28c) 

(28d) 

(28e) 

(28f) 


Af(v)=A’(v)A'(v). 

These  eigenfunctions  are  orthogonal  so  that 


(28g) 


/: 


i 


L 


dp  <Mh)<M/-0=0. 

I 

dpp  4>l±(p)=No±, 


(29a) 


(29b) 


-I 


dft-P  4>„(n)<f>r  (p)=JV(p)8(v-v'). 


(29c) 
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We  remark  here  that  this  set  of  eigenfunctions  seems  to  be 
rather  natural  for  the  expansion  off(p—*p<,)  for  the  reason 
that,  as  we  shall  see  later,  for  both  extreme  cases 
/lfi-*fto)  = l (isotropic)  and = the 

variational  principle  gives  exact  results.  Now  consider  the 
expansion, 

/V<h  ^ 

+ f (30) 

J- 1 N(v) 

so  that 

B0±(po)= J dpp/lp-^p,,)^^)  (31a) 

and 

^vG-O-J  dpp/\p-*po)<\>v(p).  (31b) 

Since/tji— *>io)  is  symmetric,  the  expansion  (30)  on  the 
right-hand  side  must  also  be  symmetric  in  p,po.  Inserting 
Eq.  (30)  in  the  integral  equation  (23)  we  get 


Ao.  J-  * W(v) 

-S»(p,po)-  4-  f dp' S„(p',po) 

2 J-i 

X -J-  [<Mm)W)~<Mp)  W)1 

+ f ^4>v0*)*v0ol-  (32) 

J-iiV(v)  1 


Upon  using  the  orthogonality  properties  (29a,  b,  c)  in  Eq. 
(32),  we  obtain 

J dp'  S,(p',p)B0±(p')  (33a) 

and 

Bv(p)=Av(p)~  -f  f dp' S,(p\p)Bv(p'),  (33b) 

2 J-i 

where 

^o±(M)=J  dp'p’ +0±(p')SJp',p)  (34a) 

J dp'p'  i()i'W#)'  (34b) 


By  comparing  integral  equations  (33a,  b)  with  Eq.  (25)  we 
see  the  obvious  parallel  and  note  that  variational  principles 
analogous  to  that  for  Eq.  ( 23 ) can  be  set  up  for  B0±  ( p ) and 
Bv(f t)  to  solve  the  inverse  problem.  As  mentioned  earlier  the 
advantage  of  dealing  with  the  integral  equations  for  expan- 
sion coefficients  B0±  and  B„  is  that  in  both  extreme  cases  of 
isotropic  scattering  l/lp-*p«)  = 1]  and  when/tp.— ►/!«)  is 
purely  forward  peaked  [monodirectional,  i.e., 

=d(fi—fio)]  the  variational  principle  provides 
exact  solutions  for/lp~+p,).  We  may  see  that  as  follows: 
Consider  Eq.  ( 33a)  for  B0±(pt)  and  let  the  trial  function  be 

♦o±(M)=«o±^o±(M)  (35) 


where  a<,±  is  a discrete  parameter.  Also  consider  the  func- 
tional for  (33a)  analogous  to  (13b) 

^±[^±1=  J ( dp  <t>o±0i)[  2/t0±(M)-4>0±(M) 

- y|‘  ^ dp'  S„(>i»<t>0±(>i')] . (36) 

Put4>0±(ft)=a0±^o±(fi) in  (36)  to  obtain 

f‘o±['do±]  =ao±(?— ao±)Bo±  ^~ao±Qo±>  (37) 

where 

B0 ±=J  dpA$±(p),  (38a) 

Qo± - f d\ t A0±(p)J'  dp' 5. (p',p)A0±(p').  ( 38b) 

Maximization  of  the  functional  F0±[^0±],  defined  by  Eq. 
(37).  yields 

-”*-(1+f^)-'  (39> 

Similarly  for  the  continuum  expansion  coefficients  Bv(  p ) of 
Eq.  (33b),  if  we  consider  the  trial  function 

<&w(p)=avAv(p),  (40) 


then  the  maximization  of  the  functional 


= J dp  <bv(M)[  2A„(p)-<PJp) 

- j- J' dp' 

(41) 

yields 

(42) 

where 

Pv=j  dpA2v{p), 

(43a) 

Gv= J’  dp  AJp)  J"  dp'  SJip',p)Av(p'). 

(43b) 

Hence,  by  replacing  B0±  by<J>0±  and  5„by  <l>„in  Eq.  (30), 
we  have  an  approximate  representation  of  the  phase  func- 
tion, which  is  given  by 

<lp-.p,) » JL  [^(Mc)  <Mn)/(i  + 

-^(Ms)^(M)/(l+Y^)] 

We  may  now  check  the  accuracy  of  the  representation 
(44)  for  the  two  extreme  cases  when 

<«> 

M.  Kanal  and  H.E.  Moms 


1796 


J.  Math.  Phys.,  Vot.  19,  No.  6.  August  1978 


1796 


rf  * 


I: 


J 


For  the  two  cases  we  conclude  from  Eq.  (25)  that 


isotropic, 

S.W,)=  2 ,-c 

8(/i— /Ao)  monodirectional. 

.1-c 

From  Eqs.  (38a,  b)  and  Eqs.  (43a,  b)  we  obtain 

Qo+  Qv  2 


for  both  isotropic  and  monodirectional  phase  functions, 
while 


4o±  — 
Av(fi)=v 


isotropic 


monodirectional. 


1—  c 


Thus  for  the  isotropic  case  we  conclude  from  Eqs.  (44)  and 
(48)  that 

=0  -c)(  ~ [<M*)+<Mm.)1  + £ _ ~j. 


But  the  right-hand  side  of  Eq.  (50)  is  merely  a unity.  This  is 
easily  seen  by  noting  that  Case’s  eigenfunctions  satisfy  the 
following  relations: 


J dfifi  1 — c), 

j*  dftfi<h0±()i)  = (±v,)(l-c). 


An  integration  of  the  completeness  relation  (27)  with  re- 
spect to  ft  gives 

<52» 

Hence  for  the  isotropic  case  we  obtain  in  Eq.  ( 50), 
/Tft-*jx,)  = 1 (as  expected).  For  the  monodirectional  phase 
function,  substitution  of  Eqs.  (49)  and  (47)  in  Eq.  (44) 
merely  reproduces  the  completeness  relation  ( 27 ) save  for  a 
factor  of  two  so  that/( ft— »fto)  = 28(  ft  — fto);  the  factor  of 
two  should  be  there  because  of  the  normalization  so  that 


ij  m>)  = i- 


In  conclusion  we  wish  to  point  out  that  for  the  direct 
problem  one  can  actually  calculate  the  coefficients  A0±(p) 
and  A„(fi)  in  the  approximate  eigenfunction  representation 
(44)  of/lfi^-ft,).  Let  us  assume  that  the  right-hand  side  of 
Eq.  (44)  is  an  exact  representation  for  some  function 
»*)  so  that 

1797  J.  Malt).  Phy*.  Vot.  19.  No.  8,  August  1979 


fa(H-*fi9)=  — (<f>o.(fi0)A  o.(n)ao.- 4>o-(fi0)A  o-(fi)ao-) 
N <w 


+ f -^T^vOxo)^v(^)a, 
J-i  N(v) 


where  a<,±,  a„  are  defined  by  (39)  and  (42)  respectively. 
This  function  then  corresponds  to  a situation  in  which  some 
(abstract)  operator  is  defined  to  have  the  spectral  density 


= 1<V<1 

v N(v)a  „ 


a<JHy— v*.)dv  a^(v— v*.)dv 


M>i. 


Thus,  one  may  use/0(  ft— »ftj,  in  analogy  with  inverse  scat- 
tering, as  a sort  of  higher  order  comparison  potential.  This 
offers  an  alternate  method  of  dealing  with  the  direct  problem 
when  the  Legendre  expansion  is  not  desirable  in  the  main 
transport  equation  when/Tft-*ft»)  is  highly  anisotropic. 

For  the  purpose  of  the  direct  problem  we  now  calculate 
the  coefficients  y40±(ft)  and  Av(p).  Thereby  we  obtain  a rela- 
tion between  the  basis  used  by  Case2  and  the  Legendre 
polynomials 


A vW=£  dfiy + v(fOs.(/t». 


From  Eq.  ( 16)  we  have 

S.(ft',ft)=  £ (M'fc.OO-  (55) 

n — 0 2 

Solving  Eq.  ( 18)  forg„(ft)  and  substituting  that  in  (55) 
yields 

S^/t',ft)=  £ (2n+ 1)  F.Ot'V’.Ot).  (56) 
»-o  1 cfH 

From  (34b)  and  (56),  we  get 

Av(p)=  £(2n+ (57) 

n-0  1 ~Cf, 

where 


K(v)=f  dpL+JWJM 


and  used  the  fact  that 

fi<h„(  ft ) = v<J>„(  ft ) — vc/2.  (59) 

Actually,  1Fb(v)  are  a special  case  of  the  orthogonal  polyno- 
mials obtained  from  the  transport  equation  by  direct  expan- 
sion of in  terms  of  Legendre  polynomials.  In  other 
words,  when /„=<5b0,  we  have  from  (58),  (59)  and  the  re- 
currence relation 

(2n  + l)ii/,B(/i)=(n  + l)/>B+10i)+n/>  *_t(ft)  (60) 

the  three-term  pure  recurrence  relation  for  IFB(v): 

(" + ,(v)+ n Wn_  ,(v)=(2i»  + 1X1 


with  lf*(v)  = 1 and  IF.  ,(v)=0.  That  Eq.  (61 ) isapure  three- 
term  recurrence  relation  implies  that  Wn(v)  are  orthogonal 
and  vice  versa.  But,  more  importantly,  If'.(v)  can  be  easily 
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calculated  either  from  Eq.  ( 6 1 ) or  from  its  explicit  represen- 
tation (58).  All  these  considerations  apply  for^0i(M)  by 
merely  replacing  v by  ± v*  where  v,  is  the  zero  of  the  disper- 
sion function  ( 28e ) . We  shall  defer  any  further  discussion  of 
the  direct  problem  until  Part  2,  which  will  be  a sequel  to  this 
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A maximal  variational  principle  is  used  to  construct  an  infinite  medium  Green's  function  for  treating  the 
boundary  value  problems  of  the  linear  transport  theory  (neutron  and  radiative).  For  the  neutrons  we 
consider  the  one-speed  case  and  correspondingly  for  the  radiative  transfer  the  monochromatic  case.  The 
scattering  properties  of  the  medium  are  presumed  to  be  dependent  on  the  relaxation  length.  Thus,  for  the 
neutrons  the  secondary  production  function  depends  on  the  neutron's  relaxation  length  and  for  the 
radiative  transfer  the  albedo  for  single  scattering  is  dependent  on  the  optical  depth.  These  two  functions 
are  kept  arbitrary  so  that  a large  class  of  problems  can  be  covered.  The  basic  principle  involves  a 
functional  which  is  an  absolute  maximum  when  the  trial  function  is  an  exact  solution  of  an  integral 
equation  of  the  Fredholm  type  The  kernel  of  the  integral  equation  is  required  to  satisfy  certain  symmetry 
and  boundedness  properties.  We  also  exhibit  an  interesting  relation  between  the  absolute  maximal  and 
Schwinger's  stationary  variational  principles,  which  in  general  is  neither  a maximum  nor  a minimum. 


1.  INTRODUCTION 

A variational  principle  for  the  solution  of  integral 
equations  with  probability  kernels  was  used  by 
Demarcus1  to  solve  the  problem  of  Knudsen  flow.  It 
involves  a functional  which  is  an  absolute  maximum 
when  the  trial  function  is  an  exact  solution  of  the 
transport  equation.  This  variational  principle  is  a 
special  case  of  one  discussed,  for  example  in  Ref.  2, 
for  which  a functional  is  an  absolute  maximum  when  the 
trial  function  is  an  exact  solution  of  an  integral  equa- 
tion of  the  Fredholm  type,  such  that  the  kernel  is  re- 
quired to  satisfy  certain  symmetry  and  boundedness 
properties.  The  more  general  variational  principle  has 
already  been  used  for  the  Gel’fand— Levitan  equa- 
tion.5,4  In  the  present  paper  we  use  the  variational 
principle  to  provide  insights  into  solutions  of  unsolved 
problems  of  transport  theory  such  as  radiative  or  neu- 
tron transfer  in  inhomogeneous  media. 

The  abstract  principle  is  discussed  for  example  in 
Mikhlin’s  book1  which  we  follow.  It  has  been  utilized  by 
Demarcus1  and  much  later  by  Stokes  and  Demarcus.’ 
Our  work  differs  from  that  of  the  others  in  the  funda- 
mental sense  that  we  use  the  principle  to  construct  the 
approximate  Green’s  function  and  to  estimate  the  spec- 
tral function  of  the  transport  operator  for  the  trial  func- 
tion chosen.  Greater  flexibility  is  thus  achieved  with 
respect  to  boundary  conditions. 

2.  THE  VARIATIONAL  PRINCIPLE 

The  basic  principle  underlying  the  variational  ap- 
proach (see  Ref.  2)  is,  that  given  the  Fredholm  integral 
equation 

f[x)=d>ix)+ £ dyK(x,y)f[y),  (1) 

the  functional 

ff«|  = /*  dx  nU)[2$(x)  + jf*  <fy  K{x,y)n(y)  - n(x)|  (2) 

is  an  absolute  maximum  when  n{x)=fi r),  for  the  kernel 
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K(x,y)  which  satisfies  the  following  conditions  of  sym- 
metry and  boundedness: 

(a)  the  kernel  K(x,y)  is  symmetric  in  x,  v, 

(b)  K(x,y)  > 0 for  x b and  a *■  y * ft, 

(c)  I dxK(x,y)  < 1 for  <i  s x «.  b. 

For  most  transport  problems  it  turns  out  that  the  in- 
tegral equation  that  one  needs  to  consider  (for  example 
for  an  infinite  medium  Green’s  function)  is  not  Eq.  (1) 
but 

fix , x0) - <t>(x , x„) 4-  I *dv’ K(x,x')fix’ ,x0),  (3) 

« 

where  v„  is  an  additional  parameter.  A further,  rather 
pleasant,  property  that  is  encountered  is  when  the  in- 
homogeneous term  <b{x,xB)  is  the  same  as  A'(v,v#).  Then 
Eq.  (3)  is  in  a symmetric  form.  That  is,  when 

<My,y0)  = K(v,x0), 

Eq.  (3)  becomes 

Av.*0>  = A'U.  »0)  + I*  dx'  Kix‘ , x)fix' , x0),  (4) 

and  the  general  functional  to  be  considered  is 
F|n|  = / * rfv  »(>  , v0)|2Mv,  v0) 

+ | * dx' K (r , x' )n(x' , v0)  - n (v , v0) | , (5) 

which  is  an  absolute  maximum  when  tt(x , x„)  - fix , y„) 
(c.f.  Refs.  2 and  3).  The  resulting  value  of  the  func- 
tional is  then 

Flfl=rdxfix,xt)Kix,x„).  (6) 

’ • 

Using  Eq.  (4)  in  |[n.  (6)  to  eliminate  the  integral  we  get 

Fl/|  = lim  [f{r,  v0)  - K (y , a0) | . (7) 

*•»« 

Thus,  for  an  exact  solution  of  the  integral  equation  (4), 
the  absolute  maximum  value  of  (he  functional  is  the  dif- 
ference of  boundary  values  of  fix,  \„)  and  the  kernel 
A"(v,  v„)  as  y approaches  from  either  side  of  v„.  Two 
further  points  are  worth  noting.  First,  in  most  tran- 
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sport  problems  of  interest  the  density  function  (fix,*, )) 
and  the  scattering  kernel  (/C(r,x0))  are  singular  at 
x = x0.  However,  the  difference  of  the  two  functions  at 
x = v0  is  bounded  from  above.  Furthermore,  as  we 
shall  see  later,  the  value  of  that  difference  is  related 
to  the  spectral  density  of  the  transport  operator  cor- 
responding to  the  choice  of  the  trial  function  [for  ex- 
ample see  Eq.  (57)  in  Sec.  3].  It  is  precisely  for  such 
a relation  and  the  maximal  property  of  the  functional 
that  it  is  possible  to  obtain  the  best  estimates  of  the 
complete  spectrum  of  the  transport  operator  and  hence 
the  Green’s  function. 

Now  we  consider  an  application  of  this  formulation 
to  the  problems  of  radiative  transfer  in  a vertically 
inhomogeneous  atmosphere  where  the  albedo  for  single 
scattering  is  an  arbitrary  function  of  the  optical  depth. 
The  same  application  is  valid,  aside  from  changing 
of  names,  for  the  neutron  transport  in  inhomogeneous 
media  if  the  corresponding  albedo  for  single  scattering 
(or  secondary  production)  of  neutrons  does  not  exceed 
unity  anywhere  in  the  given  medium,  i.e. , for  the 
subcritical  problem. 

3.  APPLICATION  TO  PROBLEMS  OF  RADIATIVE 
TRANSFER  AND  ONE-SPEED  NEUTRON  TRANSPORT 

We  shall  consider  the  problems  of  radiative  trans- 
fer. Hence  we  use  the  corresponding  appropriate  term- 
inology. The  equation  of  radiative  transfer  for  a plane 
parallel  vertically  inhomogeneous  isotropically  scat- 
tering atmosphere  is 

m 57  (x,p)+/(x,p)=  f du'Hx,n')  + Q{x,n), 

dX  2 J.X  (8) 

where  p is  the  cosine  of  the  angle  between  the  vertical 
axis  and  the  direction  of  propagation  of  radiation,  x is 
the  optical  depth,  w(x)  is  the  albedo  for  single  scat- 
tering, and  f(x,p)  is  the  specific  intensity.  We  have 
added  a source  term  Q(x,  p)  in  case  one  is  dealing  with 
neutron  transport.  For  problems  of  radiative  transfer 
Q will  be  zero  In  general.  Let  the  atmosphere  be  finite 
and  optically  bounded  between  x = xt  and  x = r, . In  order 
to  set  up  the  boundary  valve  problem , we  follow  the 
standard  procedure  discussed  in  Refs.  6 and  7.  Thus, 
consider  a time  reversed  adjoint  equation 

- M jji  GCx,  - P "*  x0,  — p0)  + G(x,  - P - x0,  — p„) 

+ ^-8(x-x0)fi(p-p0).  (9) 

Note  that,  here,  the  delta  functions  are  weighted  by 
Vw(x)/2.  Multiplying  Eq.  (8)  by  G and  Eq.  (9)  by  /, 
subtracting  and  integrating  the  difference  with  respect 
to  p from  - 1 to  + 1,  and  x from  x,  to  xt,  we  get  the 
result 

fix,  P)  = <fp'  M'[/U„  p')G(x„  - p'  - x,  - p) 

-f(x,,p')C(x,,  -p'-x,  -p)] 


x G(x* , - p'  -x,-p). 

(10) 

In  particular,  If 

VU,p)-  2 6(x-x,)6(p-  p,) 

(10') 

and  the  medium  is  infinite,  then  /(*,  p)=G(x, 

p-x„p,) 

is  a Green’s  function  for  the  plane  source  at  .v  = r,  and 
p=  P,  satisfying  Eq.  (8).  With  suitable  boundary  condi- 
tions that  both  G and  G vanish  at  x = ±°°,  we  get  from 
Eq.  (10)  the  following  reciprocity  relation, 

V( a{x)C(x , - p-x,,  - p,)=  Vulx,  )<?(*,,  p,  - x,  p). 

(10") 

One  may  use  this  reciprocity  relation  to  replace  G by 
G.  However,  we  prefer  to  use  G instead.  For  the  finite 
medium  and  any  given  incident  conditions  for  f(x,p) 
at  the  boundaries  (i.e.,  at  v = .v,,  p>0,  x=x„,  p<0), 
the  reflected  and  transmitted  radiations  are  determined 
from  the  limits  of  the  integral  identity  (10).  In  other 
words,  the  unknown  quantities  I(xt,  p)  for  p < 0 (the 
reflected  intensity)  and  /(v6,p)  for  p >0  (the  transmitted 
intensity)  may  be  determined  from  the  integral 
equations 

/U.,p)  = -^TjT  <fp'  p' [/(*,,  p')G*(.«r,, -p'-x.,-p) 


-/(*,,  P')G(x,,  - p'  -x#,  - p)] 
2 


j=^jf ** jf  W , P'JGU' , - p'  - jr. , - p) , 


p > 0 


(11) 


and 


fl*,,p)  = ) j dp'  p'[/U.,p')CCvi,  - p* -x,,  - p) 

- ftv  p’)G*(x,,  - p'  - xt,  - p)] 

+ jh?f^dX,£ 

xG(x\  - p'  — x»,  - p),  p<0,  (12) 


where  +(-.)  on  G represents  its  boundary  value  as  x 
approaches  the  boundary  of  the  atmosphere  from  the 
right  (left).  Solution  of  the  integral  equations  (11)  and 
(12)  is  possible  provided  we  can  construct  the  Green’s 
function  from  the  adjoint  Eq.  (9).  In  fact  the  proof  of 
the  existence  of  solutions  of  integral  equations  of  this 
type  for  the  grey  problem  is  guaranteed  by  construc- 
tion of  the  unknown  coefficients  (see  Ref.  6 for  further 
discussion).  We  focus  our  attention  on  procedure  for 
the  construction  of  G for  the  simple  (and  obvious)  rea- 
son that  It  is  the  central  quantity  required  for  any 
boundary  value  problems  of  this  kind.  We  shall  dem- 
onstrate here  the  use  of  the  variational  technique  to 
obtain  an  approximate  solution  to  that  problem  of  de- 
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termining  G.  To  do  that,  we  first  obtain  an  integral 
equations  for  G.  A convenient  way  is  to  take  the 
Fourier  transform  of  Eq.  (9)  with  respect  to  x.  The 
result  is 


G(x.-  P— x#>-  p0) 


= **'  e^P(-  iJttr'JwU'ljf1  dn’ 

x Glx*,-  p'-x0>-  M0)  + Vu>(*0)  6(p-  p„) 


X 


^ exp[i7?u  -Xp)l 
l-l*p 


(13) 


For  convenience  in  writing,  denote  the  angular  integral 
of  the  Green’s  function  on  the  right-hand  side  of  Eq. 

(13)  by 

i 

H^x\xa,pa)=^dpG(x,-  p~x 0,-p0).  (14) 

Then  we  have 


C(x,-  p-*x0P-  p0) 


dk 


exp(ikx) 

1 — ikp  J 


dx'  o>(x') 


x exp(-  ikx1  )H  lx'  ;x0,  P0) 


where  E,  ( I x - x'  I ) is  the  exponential  integral.  Note  that 
+ and  K are  related  by 

/«x,x0)=7|rfp*lx;x0,p).  (19) 

To  complete  the  set  of  integral  equations  that  we  need, 
we  define  the  density  p(x,x„)  by 

p{x,x0)  = f'  dpH(x;x„,p)  (20) 

and  from  (16)  we  find  that  p(v,r0)  satisfies  the  following 
integral  equation, 

v'u)  (x) plx,  x0)  = K(x,  x„)  + jH  dx'K(x', x) 

xVii t[x')p(x' ,x0),  (21) 

where  we  have  used  the  integral  relation  (19)  between 
+ and  K.  If  we  let 

Hs{x\x0,p0)  = 'lw  (x)  «(x;x0,p0)  (22) 

and 

Ps(x,x0)  = Va)(x)p(x,x0),  (23) 

then  the  integral  equations  (16)  and  (21)  assume  sym- 
metric looking  forms: 

U5 1*  »-^o>  Po)  ^X\xat  Pq) 

+ fm\dx'K(x',x)Hs  (x';x0,p0)  (24) 


V^6(p-  p0)  dk  . (15) 


and 


Ps  U . v0)  = K (x  ,x0)  + f rf.v'  Klx' , x)ps  tv' , x0) . 


(25) 


Clearly,  if  we  can  determine  Wlv;v„,  p0),  then  for  any 
given  b>(x),  G is  completely  determined. 

A.  Integral  equations  for  Hlx;x0,p0 ) and  p(x,x0 ) 

By  Integrating  Eq.  (15)  with  respect  to  p from  - 1 to 
+ 1 we  obtain 


B.  Functionals  for  Hs (x,x0,p0 ) and  pj  (x,x0 ) 

Consider  the  functionals 

= f"  dx  N(x;x0,  p0)(2+(x;x0,  p0)  - N(x;x„,  p„) 

+ L.  dx‘  KLx’  ,x)N(x’-,x0,  p0)|  (26) 


H(x;x0,  p0)  = *lx;x0,  p0) 


v^yj/x'x(x’> 


x)Vu)l x'J//(x'  ;x0,  p0), 


(16) 


where 


+(x;x0,  p„)=  ^ VwlxJwUo )J"  dk  (17) 


and 


K(jc'  ,x)=  jj-  Vuj(x')w(x)J^*  dk  expf  1*  (x  - x' ) 1J^  * ^ 

= J y/wOc'twlx)  f ~ exp(- |x-x'|/p) 

■Jo  ^ 

= JVw(x’)wU)  £,(|x -x’  |),  (18) 


and 

f M=  /’  rfv p(.r,  xa)[ ZK (x,  x0)  - »{x,  x0) 

+ /__rf.x’Alx',x)»lv',v0)).  (27) 

Since  the  kernel  K Iv’.x)  in  both  integral  equations  (24) 
and  (25)  satisfies  the  three  properties  listed  below  Eq. 
(2),  we  conclude  that  functionals  F,(.V]*  and  F(u|  are 
absolute  maxima  for  exact  solutions  of  those  integral 
equations  (see  Refs.3  and  4 for  the  proof).  Thus,  for 

N(.v,.Vq,  Pq)  ~ H § lx  :.Vq , Pq ) 

and 

«(*,  v0)  = ps(v,x0), 

the  values  of  the  corresponding  functionals  are 

= ^Kv;vn,  P0)WS lv;v0,  p„)  (28) 

and 

JlPsl  = .L  *vffl v,x0)ps(x,x0).  (29) 
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Further  use  o(  Eq.  (25)  in  Eq.  (29)  yields 

Ftp*  I = lin>  [P*  (*.  »0)  - KU,*0)).  (30) 

Equation  (30)  is  a rather  surprising  result  for  the  rea- 
son that  the  value  of  the  functional  for  the  density  bears 
such  a simple  relation  to  the  density  ps(v,.r0)  and  the 
scattering  kernel  K(x,x0).  This  is,  in  the  strict  sense, 
a relation  between  the  direct  and  the  inverse  problem. 
For  the  Gel’fand— Levitan  equation  (c.f . Ref.  3)  we  dis- 
covered a relation  identical  to  Eq.  (30)  where  it  turned 
out  that  the  result  can  be  translated  into  a theorem 
about  the  area  under  a curve  to  a given  point  x when 
considered  as  a functional  of  n(jc,xa).  This  curve  was 
given  by  the  scattering  potential  to  the  given  point  when 
the  functional  took  on  its  maximum  value.  In  that  case 
the  functional  may  thus  be  considered  as  a method  of  ob- 
taining the  scattering  potential  from  the  spectral  data 
through  the  variational  technique.  However,  for  the 
transport  problems  we  are  considering  such  a simple 
interpretation,  is  not  possible.  For  one  reason,  we  note 
that  both  ps(jr,  r0)  and  K(x,x0 ) are  singular  at  *=v0. 
However,  in  the  limit  when  jr—  v0  the  difference  of  the 
boundary  values  of  ps tv,r0)  and  K(x,x0)  is  bounded.  It 
would  be  of  value  to  determine  that  limit. 

C.  Value  of  the  functional  Ffp$  ] 

Clearly,  it  is  of  crucial  importance  to  know  the  value 
of  the  functional  F[ps ) for  the  exact  solution  of  the  in- 
tegral equation  (27).  The  reason  is  simply  that  for  any 
sequence  of  trial  functions,  the  values  of  the  corre- 
sponding functionals  must  approach  the  absolute  maxi- 
mum value  obtained  for  the  exact  solution.  One  can  then 
estimate  the  absolute  error  introduced  for  any  particular 
choice  of  a trial  function.  Of  course,  in  order  to  find 
f(ps],  as  given  by  Eq.  (30),  we  must  know  the  exact 
form  of  the  density  ps{x,x„)  (which  in  fact  is  the  un- 
known). However,  since  the  functional  is  the  boundary 
value  of  the  difference  between  ps(.r,.v0)  and  K(v, *„) 
when  x approaches  .v0,  we  may  take  advantage  of  that 
to  construct  a pseudo  Green’s  function,  such  that  if 
p<Cr,.r0)  is  the  corresponding  pseudo  density,  the  dif- 
ference ps(x,x0)-  P,(x,x0 ) in  the  limit  .v— ,r0  is  small. 

In  other  words  we  require  that  pf(x,.v0)  be  singular  in 
the  same  manner  as  ps(.v,.r„)  at  ,v  = x0  plus  a (possible) 
small  correction  term  which  is  nonsingular  at  x = x0 
for  all  values  of  x„.  We  proceed  as  follows:  Consider 
Eq.  (9)  for  the  Green’s  function 

4 G(x,  4 - x0,  40)  + G(x,  p - x0,  p,) 

= ^2^/'  *4'C(x\p'-*t,pD)  (31) 

+ fi(*-*0)fi(p-  P0)» 

where  for  convenience  we  have  changed  the  signs  of  p 
and  p0.  Now  consider  a pseudo  Green’s  function  which 
satisfies 

p|7G<(x,p-X0,Me)+Gf(x,M-r0,M0) 

= ^ J ' W C,<*.  M'  - *0,  4„)  (32) 


+ ^2^'  6k-x0)6(p-  P„). 

Note  that  u>  associated  with  the  first  term  on  the  right- 
hand  side  in  Eq.  (32)  is  a function  of  x0  (the  source 
point).  If  we  let 

L(x,  p - x0,  p0)  = G(at,  p - v0,  p0)  - Gftv,  p - x0,  p„),  (33) 

then  subtracting  Eq.  (32)  from  Eq.  (31)  we  find  that  the 
difference  L(x,  p — x0,  p„)  satisfies  the  following 
equation, 

U Jx  ^v’  M ~ *<”  ^o)  + Uv,  P - .V0,  p0) 

= Jf  dp' Ux,  p'--v0,  P») 

+ — ^ dp'  Gf(.v,  p'  - x0,p0).  (34) 

In  what  follows,  we  assume  that  u.>(.v)  is  either  a con- 
tinuous function  or  can  be  approximated  by  a suitable 
continuous  function.  If  u>  were  independent  of  v,  then 
L(x,  p — .v„,  pn)  would  be  the  infinite  medium  Green’s 
function  without  sources.  In  consequence,  with  ap- 
propriate boundary  conditions  at  x = ±^,  the  only  solu- 
tion of  Eq.  (34)  would  be  the  trivial  one.  i.e. , L =0, 
as  one  would  expect.  In  that  case  G would  coincide  with 
Gf  and  the  value  of  the  functional  fas  given  by  Eq.  (30) | 
may  be  readily  shown  to  be  proportional  to  the  spec- 
tral density.  For  an  arbitrary  u>(v) , the  last  term  in 
Eq.  (34)  provides  the  "source’’  for  L(x,  p — v„.  p0)  and, 
therefore,  that  equation  may  have  a nontrivial  solution. 
However,  we  are  only  interested  in  the  behavior  of  I. 
near  and  at  v=  v0.  The  first  simple  conclusion  is  that 
L is  continuous  everywhere  on  the  real  axis.  This  fol- 
lows from  the  jump  condition.  Thus,  if  we  integrate 
Eq.  (34)  from  left  to  right  of  v = .vn,  we  get 

p[t‘U0,  p - v0,  pn)  - L-(.v0,  p - v0,  p„)|  = 0.  (35) 

where  we  have  assumed  that  u>(v)  is  continuous  and 
L*  = liin  L(x,  p — .v0,  p0). 

*•*0  (Right  of  x0) 

(Loft  of  r0> 

Continuity  of  L also  implies  that  its  derivative  is  also 
continuous  at  v = .v„.  This  follows  from  simply  taking  the 
difference  of  boundary  values  of  Eq.  (34)  as  v — v„  from 
both  sides  of  ,v0.  If  we  solve  Eq.  (35)  for  p,  then  the 
most  general  solution  is 

*'*('„>  4 *o>  40)  - (-"('„,  p - xn,  p„)  = A06(p).  (36) 

But,  from  Eq.  (34)  for  p = 0,  we  get 

C* , 0 -•  x0 , p0)  = f dp'  L(x,  p'  — a0,  P0) 

J-\ 

+ r dp'  CfU,  p'  - -V0,  p0).  (37) 

J-l 

Thus, 

Uo i®"*  *o»  40)  * Uo i0  "*  Aot  4p)  =0i 

so  that  \f,*0  in  Eq.  (36).  We  have,  therefore,  shown 
that  L(.v,  p — .v„,  p,)  is  continuous  everywhere  and  so  is 
its  derivative. 
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We  now  proceed  as  before  as  we  did  with  Eq.  (9)  and 
obtain  the  integral  equation  for  pt(x,x0)  defined  by 

Ptl*,x0)  = /‘  duf^dMoUx,  fi~x0,p0).  I38) 

The  integral  equation  satisfied  by  pt(x,x„)  is 
PtU,*0)  = i /.  -x|)w(x')ptlx\x0) 


we  have  from  Eqs.  (30)  and  (44)  the  value  of  the 
functional 

F(pi  ] = lim  [ps  (.v,  v0)  - KCv,  *„)] 
r**o 

approximately  equal  to 

F(ps  1 a lim  [pfS  (x , x0)  - K(x , x0)  1 

"‘o 


+ {fml  rfx'E1(|x'-x|)[w(x')- a>Cv0)] 
xpf(*/,x0),  (39) 

where 

P,<J‘,x0)  = f i rf M _|T|  du0Gt{x,  p ■*  x0,  p0)  (40) 

is  the  pseudo  density  and  E,(lx'  -xl ) is  the  exponential 
integral, 

£,(|x'-x|)  = f exp(-  |x' -x|/p.  (41) 

Jo  ** 

Since  pL(x,x0)  is  continuous  at  x = x„,  we  get 
PtU0,x0)  = i dx’  £,(|x'  - x„|  )w(x')pIl(x',x0) 

= 5 /_  dx'  £>  ( | x'  - x°  | )[ w (x' ) - w (x0)  |pf  (x’ , x0) . 

(42) 

Now  we  make  the  approximation  for  the  first  integral 
on  the  right-hand  side  of  Eq.  (42).  The  largest  con- 
tribution that  the  exponential  integral  £,( lx*  - ,v0l ) makes 
is  near  x'=x„.  In  fact  £,  is  weakly  singular  (logarith- 
mically) there  and  decays  rather  rapidly  away  from  x0. 
In  analogy  with  the  method  of  stationary  phases,  we 
assume  that  pL( x',x0)  varies  slowly  near  x'  = v0  and, 
therefore,  may  be  taken  out  of  the  integral.  Equation 
(42)  then  gives  us 

Pt(*o.x0)  = H jH dx>  £ill*'  “ r«  I >l<*»(x')  “ <»tx0)| 

xP,l*’,x0)}/[l  -i  | <fx'  E,(|x*  -XolJwlx*)}. 

- (43) 


+ y/wlxj  pL(x„,r0),  (47) 

where  pt(x0,  r0)  is  given  by  Eq.  (43). 

Here  we  have  assumed  that  one  can  construct  the 
Green’s  function  G'fCv,  p — v0,  p0)  by  solving  Eq.  (32) 
and  that  the  limit  in  Eq.  (47)  exists  and  is  bounded.  We 
answer  that  question  by  merely  noting  that  the  method 
of  obtaining  solution  of  Eq.  (32)  is  in  exact  parallel  to 
the  case  where  the  medium  is  homogeneous.  In  fact, 
following  Case  and  Zweifel  [see  Eq.  (11),  p.  97,  in  Ref. 
6[,  we  find  that 


0o>(p)0<u(Po)exp(—  I > - ml  t>a) 


Gf(x,  p — .v0,  p0) 

-4  ^ 

+f  ^f^«.>^.>o)exp(-l'-  Vol  f)).  W8) 

*'o  f 

where  +(-)  is  for  x>x0  (r<  x0),  and  *o.  (p).  </>„(p)  are 
discrete  and  continuum  Case  eigenfunctions  given  by 


<t>o.(v)  = 


iwtvo) 

2(±  v0  - p)’ 


</>,( p)=  ■ -if  — p +A(p)a(i/-  p), 

A(v)  = 1 - i'ui(.v(1)tanh'1i>, 

..  U>(vn)  j/(i»(vn)  1 \ 

w T-T> 

N(v)  — i’l (1  -a>(.vJfManh‘‘r)2  + u>2(.v0)rV  4|, 


(49) 

(50) 

(51) 

(52) 

(53) 


But,  by  definition  [see  Eqs.  (33),  (38),  and  (40)| 


v<*>tx0)  Pt(x0,x0) 


= lim  ['/u>C*)p(x,x0)  - '/<*i(x0)  p,(x,x0)| 

«-IO 


therefore. 


OTptCx.,x0) 

= lim  ['/a) (x)  p(x,x0)  - )C(x,x0) 

»*«o 

- Vwlx0)pf(x,x0)  + KCx,x0)|,  (44) 

where  K(x,x0)  is,  by  definition  (18),  given  by 

E(x,x0)  = J Vw(x)w(Xo) £,(|x-x„|).  (45) 

Since,  previously  we  defined  p, (x , x„)  = Vw (x)p(x , xe) 

[see  definition  (23)]  and  now  we  define 


P„  lx,  x0)  = Vu>(x0)  p,(x , x0) , 


(46) 


and  is  the  zero  of  the  dispersion  function,  i.e.. 


U>(\o)l»0 

2 


(54) 


We  remark  here  that  Cases’s  discrete  spectrum  ob- 
tained from  the  solution  of  Eq.  (54)  is  in  fact  now  a 
“Regge”  type  trajectory  determined  by  the  functional 
form  of  w(xn).  Further,  the  functional  F(ps|  will  also 
follow  a trajectory,  which  we  will  call  “maximal  tra- 
jectory." It  is  now  a simple  matter  to  find  that  maxi- 
mal trajectory.  This  we  do  by  noting  that  tf>„,(p)  and 
$„(p)  are  normalized  to  unity,  i.e. , 

f'^dp  <t>„,(p)  = /(  dp  4>„(p)  = 1 (55) 


so  that 

, , u>Uo)/exp(-  I x-x0l  ,■><,) 

Pf*U,x0)  = ^ 

+jT  ^y«p(*'x-x,l/p)). 


(56) 
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Hence 


MpJst  ^[k  +!i?jf  dl’ex«-'x'/l’)(m  ~ :)] 

+ Vu)Cx0)  PL(x0,x0), 


where,  as  a reminder,  pL (x0,x0)  is  given  by  Eq.  (43). 
From  Eq.  (57)  we  see  that  the  functional  is  bounded.  In 
the  lowest  approximation  the  continuum  part,  involving 
the  integral,  may  be  neglected;  because,  for  small  x 
the  largest  contribution  from  the  integral  involving 
N(v)  comes  from  near  v = 0,  where  N(i>)~v.  Hence 


f°[pJ=  +Vu >(x0)ptU0,x0), 


(58) 


where  we  have  subscripted  the  functional  with  zero  to 
indicate  the  zeroth  approximation.  We  close  this  sec- 
tion with  one  more  remark.  Here  we  have  demonstrated 
the  use  of  a maximal  variational  principle  to  construct 
an  appropriate  Green’s  function  for  inhomogeneous 
media  in  dealing  with  problems  of  radiative  and  neutron 
transport.  We  have  made  no  attempt  to  solve  any  actual 
problem  for  a real  medium.  We  hope  to  make  that  a 
topic  for  the  next  paper.  However,  we  present  a simple 
example  in  Sec.  14,  where  we  consider  a trial  function 
involving  Gf. 

We  also  wish  to  point  out  that  the  same  variational 
principle  seems  to  be  applicable  to  wider  variety  of 
problems  in  theoretical  physics,  viz.,  kinetic  theory 
of  gases,  plasma  physics  (such  as  the  Vlasov  equation), 
and  electron  transport  in  the  upper  atmosphere. 

4.  EXAMPLE 

For  the  trial  Green’s  function  we  take  the  simplest 
form  which  corresponds  to  the  scaling  of  Gf,  i.e. , we 
let 


Gr(x,  P "*  Xg,  p0)  = o 


u>Cv0) 

w(x) 


G,U,  p - x0,  p^), 


(55) 


where  a is  a parameter  to  be  determined  from  maxi- 
mization of  the  functional  (27).  Now  it  follows  that  the 
trial  density  function  is 


pr(x,x8)  = a ^^-pf(*,x0). 


(60) 


It  is  seen,  by  replacing  wlr)  by  u>lv„)  everywhere  in  Eq. 
(25),  that  ptS  = ■An (x0)  pf  satisfies  the  following  integral 
equation 

P, s U , x0)  = K,  U . *0)  + jf  . dx’  K,  (*' , x)pt s (x 1 ’ , x0) , (61 ) 


where  the  kernel  Kt  is 

KfU,*.)=^  £,(|*-*0|). 


(62) 


Now  in  the  expression  (27)  for  the  density  functional  we 
set  n = pTS  = pT(x,xa)  and  compute  the  value  of  the 
functional.  One  may  readily  show  via  use  of  integral 
equation  (61 ) that  the  value  of  the  functional  is 


f [Prs  1 = n (2  ~ 0 >^[p,s  1 “ 01  » 

where  by  definition 


(63) 
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^tP*s)  = /,  dxK'byXc )ptSU,x„) 


and 


(57)  W*o >=/>($$  -1>„U.*0). 


(64) 


(65) 


Setting 


g[Prs)=° 


we  get 


(66) 


oUo)-i  + wu0)'fTp^]- 

Note  that  F[ptS  ] is  the  exact  maximum  value  of  the  func- 
tional for  the  pseudo  grey  medium.  The  value  of  the 
functional  F[prs  1 is  given  by 

FlPrsl=  F[pJ  + lix0)  ‘ . (67) 

To  test  the  accuracy  of  the  trial  function  chosen,  we 
compute 


f(y  ) = flgrs  1 r 
U°  f\pTs\ 


(68) 


and  see  how  far  tlx0)  is  away  from  zero.  Thus,  for  the 
present  example  considered  f(v0)  represents  the  local 
deviation  from  the  grey  approximation. 

5.  A RELATION  BETWEEN  THE  MAXIMAL 
AND  SCHWINGER'S  VARIATIONAL  PRINCIPLES 

In  the  previous  sections  we  have  illustrated  the  basic 
ideas  underlying  the  maximal  variational  principle  as 
applied  to  transport  theory.  Because  the  functional  used 
is  an  absolute  maximum  for  the  exact  solution  of  the 
appropriate  transport  equation,  it  has  the  advantr 
that,  in  principle  at  least,  one  can  obtain  the  est,  t 
of  absolute  errors  for  any  set  of  trial  functions  used. 

We  have  already  provided  an  approximate  uniform  bound 
to  the  value  of  the  functional  which  for  most  situations 
should  be  useful.  By  contrast  other  variational  princi- 
ples, such  as  Schwinger’s  or  Rayleigh  Ritz  principles, 
which  require  stationarity  only,  the  estimate  of  errors 
for  a set  of  trial  functions  for  the  entire  medium  (i.e. , 
the  uniform  error)  is  usually  not  feasible.  In  this  sense 
those  variational  principles  are  rather  limited  and  often 
very  cumbersome  when  applied  to  real  problems.  How- 
ever, it  would  be  of  pedagogical  value  to  show  a relation 
between  the  maximal  principle  which  we  have  used  and 
Schwinger’s  stationary  principle. 

Schwinger’s  variational  principle*  states  that,  given 
the  Fredholm  integral  equation 

y(r,r„)  = K(r,r0)+  f"rfr'K(*-',r-).v(r'),  (69) 

* 0 

where  the  kernel  K(r'  ,r)  is  symmetric,  then  the 
functional 

S[n|  ={  X drn(r,rt)[n(r,  r0)  - jf*  dr>  A'(r',r)«(r',r0)I}/ 

[ f drn(r,r0)K(r,ra )]*  (70) 
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is  stationary  when  n(r,r0)  = y(r,r0).  The  value  of  the 
functional  is  then 


feller  University  for  reading  the  manuscript  and  for 
making  helpful  suggestions. 


Styl  I0"drK(r,r0)y(r,r0)' 


Let 


/"  drK(r,r0)n(r,r0). 


(71) 


(72) 


By  comparing  Schwinger’s  functional  (70)  and  the  maxi- 
mal functional  F[m]  (5),  we  conclude  that  the  two  func- 
tionals bear  the  following  relationship  with  each  other, 

F[n]  = 2N  (r„)  — S[n  J>/*(r0) . (73) 

For  the  exact  solution  of  Eq.  (69),  we  have  already 
from  Eq.  (70)  thatS[y]=  l/N(r0).  Therefore, 


fljpJSfy  1 = 1. 


(74) 


In  other  words,  for  the  exact  solution,  the  product  of 
maximal  functional  and  Schwinger’s  stationary  functional 
is  unity. 
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ABSTRACT 


In  this  paper  we  deal  with  the  inverse  problem  in  radia- 
tive transfer,  which  is  equally  applicable  to  the  neutron 
transport,  for  a finite  homogeneous  medium.  We  give  a method 
for  computing  the  scattering  function  and  the  albedo  for 
single  scattering.  The  solution  is  given  in  terms  of  Legendre 
expansion  of  the  scattering  function.  A decomposition  of  the 
equation  of  transfer  is  also  given  in  which  the  relation  be- 


tween the  direct  and  the  inverse  problem  is  exhibited  via  the 


principles  of  invariance.  A relation  of  this  work  with  Case's 


method  of  approach  is  outlined.  This  work  should  be  of  prac- 
tical value  to  the  problems  associated  with  remote  sensing  of 
the  terrestrial  atmosphere  and  the  neutron  piles. 


1.0  INTRODUCTION 


In  an  earlier  paper  by  Kanal  and  Moses  [1],  an  inverse 
problem  in  neutron  (or  radiative)  transport  theory  was  solved 
for  an  infinite  medium  containing  a plane  source  emitting 
neutrons  (or  photons)  in  a direction  whose  cosine  of  the  angle 
was  fixed  arbitrarily  at  p = The  objective  was  to  recon- 

struct the  scattering  kernel  f(p'-^p)  (or  the  phase  function) 
from  the  information  obtained  from  the  experiment  which  mea- 
sured the  neutron  (or  photon)  column  density  for  all  angles; 
which  is  the  zeroth  moment  of  the  angular  density  (or  the 
specific  intensity  of  the  radiation  field).  In  contrast  the 
object  of  the  usual  direct  problem  is  to  find  the  distribution 
function  (angular  density  of  neutrons  or  the  specific  inten- 
sity of  radiation)  for  a given  scattering  function.  In  this 
paper  we  wish  to  treat  the  inverse-albedo  problem  of  radiative 
transfer  in  an  atmosphere  which  is  finite  in  the  total  optical 
depth  scale.  The  direct  albedo  problem  for  a finite  atmosphere 
is  defined  as  follows:  we  are  given  an  incident  radiation 
field  at  x = 0,  where  x is  the  optical  depth.  Also  given  are 
the  phase  function,  the  albedo  for  single  scattering  and  the 
reflection  boundary  condition  at  x = x^»  where  x^  is  the  total 
optical  thickness  of  the  atmosphere.  The  problem  is  to  deter- 
mine the  radiation  field  everywhere  and,  in  consequence,  the 
reflected  component  at  x = 0 and  the  transmitted  component  at 
x = x^.  In  contrast  the  inverse-albedo  problem  for  the  finite 


atmosphere  involves  the  construction  of  the  phase  function 
and  the  albedo  for  single  scattering  from  the  results  of  an 
experiment  which  measures  a certain  minimum  set  of  suitable 
quantities.  We  discuss  the  nature  of  measurements  that  are 
required  for  the  present  treatment  below. 

In  an  earlier  paper  [2],  K.  M.  Case  gives  an  elegant 
solution  of  an  inverse  problem  for  an  infinite  medium  contain- 
ing an  isotropic  plane  neutron  source  and  constructs  the  scat- 
tering function  from  the  spectral  data  of  the  transport  oper- 
ator. He  requires  measurements  of  the  neutron  density  for  all 
relaxation  lengths*  and  relates  the  density  to  the  spectral 
function.  He,  then  proceeds  to  show  by  different  routes  (in- 
cluding the  discretized  version  of  the  Gelfand-Levitan  equa- 
tion) how  one  may  construct  the  coefficients  in  the  Legendre 
polynomial  expansion  of  the  scattering  function.  In  our  treat- 
ment of  the  albedo -inverse  problem,  which  physically  differs 
from  the  infinite  medium  problem,  the  measurements  that  we  re- 
quire are  different  from  the  ones  required  by  Case.  In  paral- 
lel with  the  infinite  medium  problem  [1],  we  require  the  mea- 
surements to  give  information  on  the  finite  integral  of  the 
specific  intensity,  i.e.,  we  require  to  know 

x, 

J dx  I(x,p) 

X1 

and  I(x^,y),  I(x?,y),  where  x^  and  x2  are  two  arbitrary  points 
in  the  optical  depth  scale.  We  shall  then  see  how  that  quan- 
tity is  related  to  the  infinite  medium  problem  via  the  princi- 


•>*;•** *? 


j n 
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pies  of  invariance  [3],  We  shall  also  demonstrate  how  one  may 
actually  construct  the  spectral  data  for  the  Case»model  from 
our  experiment  and  in  consequence  generate  Case-eigenfunctions 
from  that  experiment. 

A final  remark  is  due  here.  In  our  previous  i )blem  [1], 
we  used  two  different  bases  for  the  expansion  of  the  phase 
function.  One  was  the  Legendre  polynomia  and  the  second  was 
the  complete  set  of  Case's  eigenfunctions  [4]  of  the  transport 
operator  for  the  isotropic  scattering  medium.  Both  expansions 
have  their  inherent  advantages.  In  the  Legendre  polynomial 
expansion,  the  inverse  problem  is  exactly  soluble  and  is  par- 
ticularly useful  when  the  medium  is  not  too  anisotropically 
scattering.  On  the  other  hand,  if  the  medium  is  highly  aniso- 
tropically scattering,  then  the  eigenfunction  expansion  is 
particularly  suitable.  For  in  that  case  we  have  a maximum 
variational  principle  which  can  be  used  (and  in  fact  we  have 
used  El])  to  obtain  a bounded  estimate  of  the  expansion  co- 
efficients. However,  the  primary  aim  of  this  paper  is  to  re- 
late the  infinite  medium  solution  to  the  inverse-albedo  prob- 
lem for  a finite  atmosphere  by  an  alternate  formulation  of 
the  principles  of  invariance.  We  shall  address  the  albedo 
problem  in  radiative  transfer. 
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2.0  INVERSE-ALBEDO  PROBLEM  FOR  A FINITE  ATMOSPHERE 


Consider  a finite  medium  optically  bounded  between  x = 0 
and  x = as  shown  in  Fig.  1.  For  the  one-dimensional,  ro- 
tationally  invariant  case,  the  standard  equation  of  radiative 
transfer  is  [3], 

|i  + = f / dp'  f(p'-p)  I(x,p ' ) . (2.1) 

Where  p is  the  cosine  of  the  angle  corresponding  to  the  unit 
vector  pointing  in  the  direction  of  propagation  of  radiation, 
x is  the  optical  depth,  I(x,p)  is  the  specific  intensity  of 
the  radiation  field,  c is  the  albedo  for  single  scattering 
and  f(p'-*-p)  is  the  phase  function.  The  phase  function  is 
normalized  to  unity  so  that 


i-  / dp'  f(p'-*-p)  = 1 . 
-1 


(2.2) 


Now  imagine  the  finite  medium,  under  consideration,  imbedded 
in  an  infinite  medium  with  the  same  scattering  properties. 
Further,  .imagine  a plane  source  of  radiation  at  x = Xq  in  the 
infinite  medium  emitting  photons  in  a direction  p = Pg.  The 
equation  of  transfer  for  that  case  is 
3f(x,p-*>xQ,pQ)  1 

M ^ ♦ V = j / dp’  f(p'-*-p)  V(x,p'-*-x0,p0) 


+ 6(x-xq)  6(p-pQ)  . 


(2.3) 


Let  us  consider  the  column  density  of  the  radiation  field  for 
the  infinite  medium  so  that. 
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FIGURE  1 


Also,  let  us  define  the  total  emission  term  by 


With  the  appropriate  boundary  conditions  at  x 


we  find  from  Eq.  (2.3)  that 


By  multiplying  Eq.  (2.7)  on  both  sides  with  f(y-^u0)  and  inte- 
grating with  respect  to  u,  we  find  that  the  total  emission 
term  Sq(|j,|Jq)  satisfies  the  integral  equation  of  the  Fredholm 
type: 


S (u,p  ) = f(p-m  ) + 


This  is  an  integral  equation  for  SQ(ii,p0)  for  the  infinite 
medium  whose  solutions  are  discussed  in  ref.  [ID.  In  the  same 


manner  as  above,  consider  the  measurement  of  the  column  density 


for  the  finite  medium  between  any  two  points  x 


where 


1 

SfCw,w1|x1,x2)  = / dp'  f(u'-ni)  Mf(p',V1|x1,x2) 


(2.11) 


is  the  emission  term  for  a slab  of  the  atmosphere  between 
x = x^  and  x = x2«  The  question  we  ask  is,  what  is  the  rela- 
tion between  the  total  emission  term  Sq(p,Pq)  and  the  emission 
term  S^(p ,p^ | x^ ,x2 ) for  the  slab,  if  we  know  the  specific  in- 
tensity I(x,p)  of  radiation  at  the  two  points  x = x^  and 
x = x2?  From  that  knowledge,  how  do  we  construe  the  scatter- 
ing properties  of  the  finite  medium?  Answers  to  these  ques- 
tions are  readily  provided  by  the  principles  of  invariance 
[3],  stated  as  follows:  "The  law  of  scattering  by  a finite 
homogeneous  atmosphere  must  be  invariant  to  the  addition  (or 
subtraction)  of  layers  of  arbitrary  thickness  to  (or  from)  the 
homogeneous  atmosphere." 

Let  us  see  how  we  may  apply  this  principle  to  solve  the 
inverse  albedo  problem.  Consider  Eq.  (2.10)  and  multiply  it 
by  f(p-*-p2)  on  both  sides.  An  integration  with  respect  to  p, 
then  yields 

1 

Sf(p,P1|x1,x2>  = / dp'p*  f(p,-►p)[I(x1,p, ) - I(x2,p’)3 

1 

♦ J / dp'  f(p'-ni)  Sf(p' ,p1|x1,x2)  (2.12) 

where  we  have  renamed  the  variables.  As  an  intermediate 
step,  expand  f(p'-*p)  in  terms  of  Legendre  polynomia. 


09 

f(p*-p)  = l (2n+l)  f PnCp')  Pn(p>  . 
n=0 


(2.13) 


I [j 


Insert  (2.13)  in  (2.12)  and  use  the  orthogonality  property  of 


Pn(p)  to  obtain. 


/ dM  pn(P>  Sf (m,m1|x1,x2) 


n " 1 


-.(2.14) 


dp  Pn(p)[2p(I(x1,p)-I(x2,M))*c  Sf  (p,px  Ix^Xj)] 

Now  in  ref.  t 1 ] (see  Eq.  (18)  in  that  reference),  we  obtained 


the  following  relation  for  an  infinite  atmosphere: 

, . W 


where 


so  that 


n " 2Pn(p0)+c  qn(yQ)  * 


qn(M0)  = [x  dw’  Fn(w)  S0(M»lJ0)  » 


s l ^ W Pn(M) 
n=0 


(2.15) 


(2.16) 


(2.17) 


Now  we  apply  the  principles  of  invariance  by  equating  the 
right-hand-sides  of  Eqs.  (2.14)  and  (2.15).  Transposing  the 
denominator  terms  to  the  corresponding  sides  we  obtain 
1 

I dp  Pn(p)  Sf (PjPj |xi,x2>  [2Pn(pQ)  + c qn<M0>] 


= / dp  qn(M0)  Pn(p)[2p(I(x1,p)-I(x2,p))  ♦ c Sf  (p.p-Jx^Xj)] 

(2.18) 

Note  that  the  second  term  in  the  bracket  on  this  left-hand- 
side  cancels  the  last  term  in  the  bracket  on  the  right-hand- 
side.  Now  multiply  Eq.  (2.18)  by  (2n+l)/2  and  sum  over  n. 


- 
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Noting  the  completeness  relation  for  the  Legendre  polynomia 


Jo  nr1  w - *0>  • 

and  the  definition  (2.17)  of  Sg(p,pg),  we  get 

1 

Sf<u0*pllxl*x2)  = f 1 dpp  S0<U,M0)  ^ 1 ^ x]_ » M ) - I(x2 


(2.19) 


,p)]  .(2.20) 


Equation  (2.20)  represents  a mathematical  form  of  the  princi- 
ples of  invariance  relating  the  emission  term  for  a finite 
slab  embedded  in  an  infinite  atmosphere  to  its  total  emission 
term  S0(p,pQ).  Thus  SQ  plays  the  role  of  a kind  of  Green's 
function.  The  symbol  \i^  is  reserved  for  the  direction  of  the 
incident  beam  at  x = 0 for  the  finite  atmosphere. 

From  Eqs.  (2.14)  and  (2.10)  we  readily  see  that  the  ex- 
pansion coefficients  f^  can  be  obtained  from  the  measurement 
M£  (p  ,p£  | x^  ,x2  ) of  the  column  density  between  x^Xg.  Thus,  we 
have 


_ 1 


n 


/ dp  Pn(p)  p ( I (x£ ,p ) - I(x2,p))'i 

1 - 1 
/ dp  Pn(p)  Mf (p,p1|x1,x2) 


-1 


For  later  pruposes  we  introduce  the  quantity 

g — 1 — cf 
n n 


where. 


1 

/ dp  Pn(p)  pdCx^jp)  - I(x2,p)) 

gn  - - 

/ dM  pn(M>  M£(p,p1|x1,X2) 


(2.21) 


(2.22) 


(2.23) 
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For  n = 0,  = 1,  we  obtain  the  value  of  the  albedo  for 

single  scattering  c. 


c = 1 - 


/ dy  y(I(x1,y)  - I(x2,y)) 
_ 

J ^ dy  Mf  (yjyj^  |x1,x2> 


(2.24) 


Principles  of  invariance,  represented  by  Eq.  (2.20)  gives 
rise  to  an  interesting  decomposition  of  the  equation  of  trans- 
fer (2.1).  For  in  the  limit  x^  -►  x2  (x2  = x),  Eq.  (2.20) 


becomes 


Sf(|i  ,y Jx-.x)  i 3I(x  u) 

Lim  - x - x = - / V^O* 

X^X  * Xi  -1  U U dX 


(2.25) 


while  from  Eqs.  (2.10)  and  (2.9)  we  obtain. 


c [I(x»w0>  + y0 


3I(x,yQ)  Sf (M0*yllxl,x) 

3 = Lim  — 


3x 


x^+x 


x - x. 


(2.26) 


Hence  from  (2.25)  and  (2.26)  we  get 

M ~ + I(x,y)  = - j / dy*  Sn(y*,y')  y*  .(2.27) 


-1 


Equating  the  right-hand-side  of  this  equation  to  the  right- 
hand-side  of  Eq.  (2.1)  we  also  obtain  the  relation 

/ dy*  Sn(y*,y)  y’  + / dy*  f(y'-y)  I(x,y*)  = 0 . 

(2.28) 

It  is  clear  that  Eqs.  (2.27)  and  (2.28)  are  completely  equiva- 
lent to  the  original  equation  (2.1)  of  radiative  transfer. 
However,  we  remark  that  Eq.  (2.28)  is  strictly  a relation  be- 
tween the  direct  and  the  inverse  problems.  Of  course,  one 


would  desire  a converse  relation  to  Eq.  (2.20)  in  which  the 
infinite  medium  measurement  is  related  to  the  finite  medium. 
Such  a relation  is  readily  obtained  from  Eq.  (2.18),  which 


gives  the  expansion  coefficients  qn(pn)  of  Sn(p,pn)  (see  Eq. 
(2.17)).  Thus, 


Pn(pQ)  / dp  Pn(p)  Sf (p,p1|x1,x2) 
= I — 

/ dp  Pfi(p)  p(I(x1,p)  - I(x2,p)) 


(2.29) 


or  in  terms  of  M^, 


W = I W 


/ dp  Pn(P>  Mf (p,p1|x1,x2) 


- 1 


7 dp  Pn(p>  p(I(x1,p)  - I(x2,p))  (2.30) 


Comparing  (2.30)  with  (2.21)  we  find  that 

2f 


n 


W = W 


(2.31) 


n 


From  this  and  Eq.  (2.17)  it  clearly  follows  that 


S0(u,M0>  = I (2n+l)  jJj-  Pn(u>  P„(M0)  <2.32) 

where  f s are  obtained  from  the  measurements  as  per  Eq.  (2.21). 
n 

Thus,  relations  (2.20)  and  (2.32)  represent  a reciprocity  be- 
tween the  measurements  of  a finite  medium  relating  to  the  mea- 
surements of  the  infinite  medium. 

Next  we  shall  consider  the  relationship  between  our 
process  of  measurement  to  the  one  considered  by  Case  [23  and 
actually  see  how  one  may  generate  Case-eigenfunctions  from  the 


3.0  RELATION  OF  INVERSE-ALBEDO  PROBLEM  TO  CASE'S  MODEL 


u 

u 

0 

[1 

0 

! 

0 

f] 

u 

0 

0 

[] 


For  the  sake  of  clarification  some  repetition  of  Case's 
work  is  essential  here.  In  his  model  one  has  an  infinite 
medium  with  a plane  neutron  source  emitting  neutrons  isotropi- 
cally at  x = 0.  The  homogeneous  equation  of  transfer  is  then. 


It  + 


ax 


+ = C 


CO  X 

l ^iif,  Pf(y)  / du'  P,(y»)  *(11') 
1=0  1 1 1 -1 


(3. 1C) 


where  the  scattering  kernel  has  been  expanded  in  terms  of 
Legendre  polynomia  and  we  have  suffixed  the  equations  with  C 
to  acknowledge  Case's  work.  Looking  for  infinite  medium  solu- 
tions of  this  equation  in  the  form  $y(y)  exp(-x/v)  one  finds 

<v-u)  *v(u)  = |^M(u,v)  (3.20 


where 


M(p,v)  = l (21+1)  f P (y)  h (v) 
1=0 


(3.30 


with 


h,(v)  = J dy  +v(y)  P£(y) 


-1 


(3.40 


and 


are 


and 


h0  * 1 • 

It  has  been  shown  that  a complete  set  of  eigenfunctions 


♦v(y)  = j cv  P[M(y ,v)/ (v-y)]  ♦ A6(v-y)  (3.50 
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A(v^)  = 0,  and  + (-)  represents  the  boundary  value  of  A as  the 
branch  cut  is  approached  from  the  top  (bottom).  Under  the 
given  assumptions,  there  are  only  a finite  number  of  real, 
simple  v.  with  |v. | >1.  They  occur  in  equal  but  opposite 


(u)  = N(v)  6(v-v') 


/ dyy  4.(y)  $.(y)  = N.  6(i,j)  (3. 10C 

_1  1 J 1 

As  an  application  of  the  orthogonality  and  completeness 


properties  one  can  calculate  total  density  corresponding  to  a 
unit  plane  source  at  x = 0.  This  is  the  solution  of  Eq.  (3.10 
when  the  inhomogeneous  term  6(x)/2ir  is  added  on  the  right. 


Then  Case  states  that,  given  measurements  of  *(x),  one  knows 
quite  a bit  about  the  spectral  function  p(v). 

The  functions  hf(v)  obey  a three  terms  pure  recurrence 


(1+1)  h 


(v)  + Jlh,  , (v)  = (21+1)  (1-cf , ) vh.(v)  = 0 . (3. 14C) 


Case's  version  of  the  inverse  problem  is  thus  reduced  to 
the  following:  Given  the  information  about  dp(v),  which  Eq. 

(3. 11C)  provides,  what  is  g(n)  = 1 - cfn? 

Now  Case  has  shown  [4]  that  the  normalization  coefficients 


persion  function  ACv)  by 


respectively.  Clearly  therefore,  the  knowledge  of  the  dis- 
persion function  is  all  one  needs  to  generate  N(v),  and  in 
consequence  the  spectral  function,  given  by  Eqs.  (3.12C)  and 
(3.130  and  all  eigenfunctions.  From  Eq.  (3.110  it  is  not 
obvious  as  to  how  one  may  in  practice  actually  obtain  the  in- 
formation about  the  spectral  function  from  the  knowledge  of 
t(x).  However,  in  principle  that  is  possible.  Our  main  in- 
terest here  is  that  from  the  solution  of  the  inverse-albedo 


problem  we  have  presented  in  the  earlier  section,  the  values 
of  f^  are  readily  obtained  from  Eq.  (2.21).  Then,  it  is  quite 
clear  that  all  quantities  in  Case's  model  which  involve  f, 's 


are  explicitly  calculable.  In  particular,  the  inverse-albedo 
problem  provides  a method  for  constructing  Case's  eigenfunc- 
tions from  an  experiment.  Further,  it  provides  an  experimental 
test  of  the  discretized  version  of  the  Gelfand-Levitan  [5] 
equation  and  all  equations  given  by  Case,  which  involve  f^, 
become  identities. 

It  would  be  of  pedagogical  valve  to  find  a relation  be- 
tween the  direct-inverse  problems  which  is  analogous  to  Eq. 
(2.28)  and  involves  Case's  eigenfunctions. 

mmV/  \) 

We  make  the  ansatz  in  Eq.  (2.28),  I(x,y)  = e $v(y) 
and  obtain 

1 1 
/ dp'  SQ(y ' ,p)  y'4»v(w')  = v / dp'  f(y'-ni) 

(3.17) 

for  all  values  of  v.  But  these  are  nothing  but  the  expansion 
coefficients  of  S0(y',y)  and  f(y'-^y)  in  Case's  eigenfunction 
base.  More  importantly,  however,  the  invariance  principles 
give  this  particular  relation  from  which  one  can  solve  for 
f(y'-*-y),  for  a given  S0(y',u)  and  vice  versa  by  use  of^he 
completeness  relation  of  the  eigenfunctions. 

It  appears  to  us  that  Case's  eigenfunctions  may  have  a 
natural  interpretation.  In  other  words,  could  these  eigen- 
functions correspond  to  some  pure  state  such  that  given  a pure 
normal  mode  of  radiation  (or  neutrons),  is  there  a scattering 
law  which  conserves  this  mode  without  dispersion?  If  so  then 
one  could  construct  a medium  which  will  pass  only  one  normal 
mode  and  not  the  others.  In  a sense  that  would  be  a filter. 


We  close  this  paper  with  a remark  that  the  inverse-albedo 
problem  for  a finite  homogeneous  medium  provides  a very  simple 
method  of  determining  the  scattering  properties  of  the  medium, 
which  appears  to  be  practical  for  both  laboratory  and  space 
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